Combined local-density and dynamical mean field theory calculations 
for the compressed lanthanides Ce, Pr, and Nd 
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This paper reports calculations for compressed Ce (4/ 1 ), Pr (if 2 ), and Nd (4/ 3 ) using a 
combination of the local-density approximation (LDA) and dynamical mean field theory (DMFT), 
or LDA+DMFT. The 4/ moment, spectra, and the total energy among other properties are 
examined as functions of volume and atomic number for an assumed face-centered cubic (fee) 
structure. These materials are seen to be strongly localized at ambient pressure and for compres- 
sions up through the experimentally observed fee phases (7 phase for Ce), in the sense of having 
fully formed Hund's rules moments and little 4/ spectral weight at the Fermi level. Subsequent 
compression for all three lanthanides brings about significant deviation of the moments from their 
Hund's rules values, a growing Kondo resonance at the Fermi level, an associated softening in the 
total energy, and quenching of the spin orbit since the Kondo resonance is of mixed spin-orbit 
character while the lower Hubbard band is predominantly j = 5/2. While the most dramatic 
changes for Ce occur within the two-phase region of the 7-a volume collapse transition, as found 
in earlier work, those for Pr and Nd occur within the volume range of the experimentally observed 
distorted fee (dfee) phase, which is therefore seen here as transitional and not part of the localized 
trivalent lanthanide sequence. The experimentally observed collapse to the a-U structure in Pr 
occurs only on further compression, and no such collapse is found in Nd. These lanthanides start 
closer to the localized limit for increasing atomic number, and so the theoretical signatures noted 
above are also offset to smaller volume as well, which is possibly related to the measured sys- 
tcmatics of the size of the volume collapse being 15%, 9%, and none for Ce, Pr, and Nd, respectively. 
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PACS numbers: 71.27.+a, 71.20.Eh, 75.20.Hr 



I. INTRODUCTION 



A number of the trivalent lanthanides undergo first 
order phase transitions under pressure that are charac- 
terized by unusually large volume changes, Ce (15%) ) 1 ' 2 ' 3 
Pr (9%)fMZ Gd (5%), 8 and Dy (6%) A These "volume 
collapse" transitions are believed to be caused by changes 
in the degree of 4/ electron correlation from strongly cor- 
related (localized) at pressures below the transitions to 
more weakly correlated (itinerant) above ,i2iiAii2i Associ- 
ated with these two regimes are differences in physical 
properties, with high-symmetry crystal structures below 
the collapse transitions characteristic of metals without 
/ electrons, and low-symmetry early-actinide-like struc- 
tures at pressures above, suggesting a greater participa- 
tion of / electrons in the bonding. The magnetic sus- 
ceptibility is of Curie- Weiss character in the localized 
regime characteristic of Hund's rules 4/ moments, and 
is believed to become temperature independent similar 
to Pauli paramagnetism in the itinerant regime. The 
latter behavior has actually been observed only for the 
collapsed a phase of Ce, but is inferred for the other 
lanthanides based on analogies to the actinides. 

In contrast to these apparently general trends, it has 
been recently shown that Nd reaches the characteristic 
itinerant cv-U structure in a relatively continuous fash- 
ion without any phase transitions of unusually large vol- 
ume change^ while in an analogous actinide system, Am 
is noted to undergo two collapse transitions^ It would 
therefore appear that the evolution from localized to itin- 
erant character in these / electron metals occurs more 



generally in a continuous manner which may or may not 
be accelerated by one or more large-volume-change col- 
lapse transitions. 

It was demonstrated many years ago that the vol- 
ume collapse transitions in the lanthanides and ac- 
tinides do not involve the promotion of electrons from 
/ to other states, and it was argued that these met- 
als were instead undergoing Mott transitions in their re- 
spective 4/ and 5/ electron systems Local density 
approximate (LDA) calculations modified to enable or- 
bital polarizationjiLi&iii to include the self-interaction 
correctionji2i22i2i and to explicitly incorporate the Hub- 
bard repulsion E//; 22 ! 23 appear to provide such transitions 
for Ce and Pr. They describe large- volume magnetically- 
ordered solutions with Hubbard split 4/ bands which be- 
come unstable under compression to paramagnetic so- 
lutions where these bands are grouped near the Fermi 
level. It has been suggested, however, that the rather 
abrupt delocalization of these modified-LDA transitions 
is an artifact of the static-mean-field nature of these tech- 
niques, and that a correlated solution would allow a con- 
tinuous transfer of / spectral weight between the Hub- 
bard side bands and the Fermi level over an extended 
volume ranged This is the case for the Kondo-volume 
collapse model of the Ce transition, based on a many- 
body solution of the Anderson impurity model, which 
identifies Kondo-like screening of the local moment as 
the driving force for the transition.i2SiS2iSi More real- 
istic calculations which combine LDA with correlated 
Dynamical Mean Field Theory (DMFT)^S so called 
LDA+DMFT; 30 ! 31 ! 32 generally corroborate these charac- 
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teristics for Ce* 3 ^* 3 * 35,36 and the same technique has been 
applied to the actinides^ 

Beyond the specific case of Ce, a more general under- 
standing of the lanthanide volume-collapse transitions or 
their absence may come from putting them in context of 
the extended evolution from localized to itinerant char- 
acter which occurs as these materials are compressed. 
As a step in this direction, the present paper reports 
LDA+DMFT calculations over a wide range of volume 
for Ce (4/ 1 ), Pr (4/ 2 ), and Nd (4/ 3 ), the last as men- 
tioned being notable for the absence of a volume collapse 
transition. These calculations have been carried out as in 
earlier work on Ce^S^ except that the spin-orbit inter- 
action is now incorporated. They also assume a face cen- 
tered cubic (fee) structure throughout, although LDA es- 
timates of the structural energy differences suggest these 
are significantly smaller than the relevant contributions 
to the correlation energy, which may then be discussed 
in terms of their volume dependence as more important 
leading order effects. 

The trivalent lanthanides follow a general structural se- 
quence under pressure, namely hep — > Sm-type — > dhep 
— > fee — » dfee, where these abbreviations are for hexago- 
nal close packed (hep), double hep (dhep), and distorted 
fee (dfcc)iifi*ii*iS These are all just different stacking vari- 
ants within the family of close packed structures, except 
for dfee, a soft L-point phonon distortion of fee, 38 a phase 
which only Ce does not assume. The heaviest lanthanides 
traverse the full sequence, whereas the lighter ones begin 
part way, e.g., Ce, Pr, and Nd are all dhep at ambi- 
ent conditions, although fee Ce is also metastable there. 
Since the same sequence is seen in Y which has no nearby 
/ states^ and is theoretically understood to depend on 
the 5d occupancy^ it would appear to have no relation 
whatsoever with the / electrons, which are then pre- 
sumed to be fully localized. Indeed, in the present work 
we find evidence for fully formed Hund's rules moments 
and little 4/ spectral weight at the Fermi level through- 
out the experimentally observed stability field of the fee 
phases (7 for Ce), and for lower pressures. 

Further compression takes Ce through a 15% volume 
collapse from the 7 (fee) to the a (also fee) phase, while 
Pr undergoes a 9% collapse from the dfee to an a-U struc- 
ture, and Nd passes from the dfee through two other low- 
symmetry structures before also reaching a-U, however, 
without any large volume changes. The present work 
finds a variety of signatures to accompany this compres- 
sion regime: substantial deviation of the moments away 
from their Hund's rules values, rapid growth in 4/ spec- 
tral weight at the Fermi level (the Kondo resonance) at 
the expense of the Hubbard side bands, an associated 
softness in the total energy, and quenching of the spin 
orbit due to the mixed- j character of the Kondo reso- 
nance versus the predominant j = 5/2 lower Hubbard 
band. The most dramatic changes in these signatures 
for the case of Ce coincide with the two-phase region in 
the 7-a transition. This behavior was observed in ear- 
lier work on Ce 33 - 34 - 35 and is consistent with the Kondo 



volume collapse scenario 

In regard to signatures in the total energy, a unique 
feature of the Kondo volume collapse model is the asser- 
tion that the 7-a transition is driven by a rapid drop in 
entropy from the 7 to the a side, reflecting screening of 
the 4/ moments. The transition therefore need not oc- 
cur at T = 0j&21 as has been documented for Ce alloys^i 
implying a featureless T = total energy. Nevertheless, 
due to the interrelationships between the different ther- 
modynamic functions, and even though it is the free en- 
ergy which actually determines the phase transition, the 
finite-T total energy should show an associated softening 
in the vicinity of the phase transition, and is therefore a 
useful diagnostic^ 

The present work finds the same signatures for com- 
pressed Pr and Nd, however, their most dramatic change 
takes place over the volume range where the dfee phase is 
experimentally observed, suggesting this phase is of tran- 
sitional character and not the end member of the local- 
ized trivalent lanthanide sequence as has been assumed. 
Only on further compression is Pr observed to collapse 
into the a-U structure, and as noted Nd has no collapse. 
An important difference among the lanthanides is that 
they become more localized with increasing atomic num- 
ber, due to the increased but incompletely screened nu- 
clear charge. This is evident in the present work by the 
offset of the concurrent correlation-related signatures to 
smaller volumes from Ce to Pr and then to Nd. One con- 
sequence is that the softening effect in the energy asso- 
ciated with the growing Kondo resonance must compete 
with the remaining, dominant part of the energy which 
has a curvature that grows ever larger with decreasing 
volume. While a region of negative bulk modulus is re- 
quired for an isostructural transition as in Ce, it may be 
more generally that a low bulk modulus favors larger vol- 
ume changes in structural phase transitions, which would 
then be consistent with the decreasing size of the collapse 
from 15% (Ce), to 9% (Pr), to none (Nd). The observed 
collapse transitions in Gd and Dy would then appear in- 
consistent; however, these involve filling of the j = 7/2 
subshell which complicates matters. 

In the remainder of this paper, the theoretical methods 
are reviewed in Sec. ^ numerical results for compressed 
Ce, Pr, and Nd are presented in Sec. 1 1 1 1 1 and then a 
summary and discussion is given in Sec. IIVI 



II. THEORETICAL METHODS 

The results in this paper unless otherwise indi- 
cated have been obtained by the LDA+DMFT(QMC) 
method, which refers to the merger of the local den- 
sity approximation (LDA) with dynamical mean field 
theory (DMFT)(2Si2S in order to create a composite 
tcchnique 3 ^* 3 ^ 3 ^ which rigorously treats onsite electron 
correlations and yet retains the material realism and 
specificity of the LDA. An essential component of the 
DMFT method is solution of an auxiliary impurity prob- 
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lem which is achieved here by a quantum Monte Carlo 
(QMC) algorithm4244 The present calculations have 
been carried out as described in previous work for Ce ;? 4 ^ 
except that the spin-orbit interaction has now been added 
and a broader set of materials considered. In the remain- 
der of this section we give a brief review of the method, 
various computational details, and calculation of the mo- 
ment. 



A. The LDA+DMFT method 

The LDA contribution to the LDA+DMFT method 
is to provide an effective Hamitonian which includes all 
valence electron degrees of freedom, 

+ \ U f nifjm hifj'm'- (1) 

Here, k are Brillouin zone vectors, i are lattice sites, 
I is the orbital angular momentum, j is the total an- 
gular momentum (£±3 except just h for 1 = 0), m = 

-j, -j+V •• ,3, nifjm = c\ fjm c. lfjm , and the prime sig- 
nifies jm 7^ j'm' . The relevant valence states for the 
present lanthanide case are 6s, 6p, 5d, and 4/, and so 
the matrices £fL DA (k) are 32 x 32. 

As described elsewhere^ the H® BA (k) are orthogo- 
nalized one-electron Hamiltonian matrices obtained from 
converged linear muffin-tin orbital LDA calculations) 4 ^' 4 ^ 
in which the 4/ site energies are shifted so as to avoid 
double counting the /-/ Coulomb interaction Uf which 
is explicitly incorporated into Eq. Q). The spin-orbit 
interaction was included as a perturbation, 4 - with the 
spin-orbit coupling parameters kept fixed across the 
bands, although evaluated at the most important ener- 
gies, namely at the respective centers of gravity of the 
occupied state density for each orbital type. The resul- 
tant i?L DA (k) matrices were calculated over a grid of 
volumes for Ce, Pr, and Nd in an assumed face centered 
cubic (foe) structure, and companion LDA constrained- 
occupation calculations used to get the volume and ma- 
terial dependent Coulomb interactions Ufm& 

Dynamical mean field theory assumes a local or k- 
independent self-energy E(ia;), which together with the 
one-body part of Eq. Q provides the lattice Green func- 
tion, 

G k {tLu) = [iuj + (i — fl£ DA (k) - E(iw)] _1 . (2) 

This is solved in tandem with an auxiliary impurity prob- 
lem defined by a similar Dyson-like equation with the 
same local self energy. 

G imp M = [g(iiu)- 1 - Z(iuj)}- 1 . (3) 

The idea is to guess an initial E = E*- 1 ), calculate the 
lattice Green function Gt from Eq. identify its k 



average with G lmp , and then find the noninteracting or 
bath Green function Q = for the impurity problem 
from Eq. . The impurity problem is completely defined 
by Q, which describes its one-body part, together with 
the Coulomb interaction for its single site as in Eq. (JTJ 
without the i index. It may be solved exactly to within 
statistical uncertainties by QMC techniques. 4 ^ This gives 
a new G lmp which together with the input C/W defines 
a new self-energy E = E^ 2 ' again via Eq. ©. The cy- 
cle is then repeated by starting anew with E = E^ 2 ) in 
the lattice Green function Eq. and so on until self 
consistency. 

Quantities of interest such as the number rij of j = 5/2 
or 7/2 4/ electrons may be obtained from the lattice 
Green function 

ni^E^/ifGk^n)]^ *, (4) 

nk 

where the sums are over k vector and Matsubara fre- 
quency oj n = (2n+l)7rT, and the trace is over those 2j'+l 
4/ states of type j. The total number of 4/ electrons is 
then nf = 7^5/2 + nj/2- Similarly the energy per site for 
the effective Hamiltonian in Eq. is 

S DM FT = ^X] Tr [fiuuWGkM e^ 0+ +Ufd, 

nk 

(5) 

which involves the double occupation 
I ' 

d=- ( 6 ) 

as discussed elsewhere^ 4 *^ In principle this is an aver- 
age over the sites i, although in practice we obtain this 
quantity form the QMC impurity problem. 

To evaluate the total LDA+DMFT energy E tot includ- 
ing all core and outer electrons, we add a correction to 
the paramagnetic all-electron LDA energy -Elda 

-Etot = EldA + EdmFT — E m LDA , (7) 

which consists of the DMFT energy -Edmft from Eq. (J5J 
less an LDA-likc solution of the model Hamiltonian 
Eq. QJ, thus "model LDA" or E m LDA- To be more accu- 
rate, our LDA total energy -Elda is scalar relativisitc and 
so omits spin orbit, and therefore the correction E m LDA 
should be and is obtained from a modification of Eq. 
with spin orbit removed. Using this modified effective 
Hamiltonian, E m LDA is determined by a self-consistent 
solution of Eqs. (J2J) and (@J for n/ = n 5 / 2 + n 7 / 2 taking a 
self-energy E m LDA = Uf(nf — i). From this, the kinetic 
energy is calculated by the first term of Eq. JSJ) and the 
potential energy by ^UfUf(nf — I). 

The simple form of Eq. Q is a reminder of the 
still fairly new and evolving nature of the LDA+DMFT 
method. There are recent discussions of the formal con- 
text of the method^I as well as some sense that a GW 
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platform in place of the LDA may allow for more natural 
integration with DMFT. 48-49 Mutual self-consistency be- 
tween the LDA and DMFT parts has also been stressed, 
i.e. , feeding a different DMFT value of n/ back into a new 
LDA calculationi^Lii This is impractical in the present 
case given the expense of the QMC solutions, although 
we do provide a perturbative change in nj in going from 
the LDA to the DMFT, and this may well be sufficient. 
A step towards answering many of these questions is sim- 
ply to provide more quantitative tests against experiment 
of LDA+DMFT implementations with clearly stated ap- 
proximations, and this is the spirit of the present work. 

B. Computational details 

While the matrices Gk, H® DA (\t), and E in Eq. @ are 
all 32 x 32, it is customary to take E non-zero only within 
the interacting or 14 x 14 /-/ block, which permits the 
same reduction in the impurity problem itself. Even so, 
one must still determine seven distinct functions of lu to 
fill out E(iw) even for the full cubic point group^ and 
more otherwise. On the other hand, if we ignore crystal 
field, then the /-/ self-energy becomes 

Ejmj' m <(iw) = T,j(iu))Sjj'S mm / , (8) 

and we need only find functions for the two spin-orbit 
states j = 5/2 and 7/2. This is an excellent approxima- 
tion over most of the volume range studied here, includ- 
ing the volume-collapse regions, where the crystal field 
splitting ranges from about 1 to 10% of the spin-orbit 
splitting. At the very smallest volumes considered here, 
this ratio approaches 50%, however, it is to be empha- 
sized this omission is only in the self-energy, and all such 
effects arc reflected in H® DA (k) where hybridization dom- 
inates anyway. In practice we input a bath Green func- 
tion of similar structure, GjSjjiS mm i , into the QMC and 
then improve statistical uncertainties by averaging over 
the six j = 5/2 and eight 7/2 states in determining the 
output G" np . 

It was noted previously that the DMFT(QMC) itera- 
tions could be greatly accelerated by separately converg- 
ing the leading, constant, Hartree-Fock-like part of the 
self-energy outside of the QMC part of the cycled For 
the present case this constant part is 

Sf = ln f -n j /(2j+l)}U f . (9) 

Therefore after each QMC cycle we extract just the fre- 
quency dependent part of the self-energy 

AE,M = zf MC (iu) - sf ({nf MC }) , (10) 

using here the QMC rij values 

ra? MC = (2j + l)[l + Gf p (0+)], (11) 

where G™ p (r) is the imaginary time Fourier transform of 
G™ p (zcj). Then with AT,j(iuj) held fixed, we construct a 



new self-energy T,j(ioj) = E^- ({rij}) + AEj(iu>) but this 
time taking the rij defining E^ ) as obtained from the 
lattice Green function via Eq. J3J. Inserting this self- 
energy back into the lattice Green function Eq. leads 
to a self-consistent condition on both Gk(iw) and the 
rij which may be completely converged with negligible 
cost, again with the AEj held fixed. Only after this con- 
vergence would we take the resultant Gk(iuj) and begin 
the process of mapping back onto the auxiliary impurity 
problem in order to update AEj (iui) . 

We found the procedure of converging E~. in between 
each QMC iteration for AE(iw) to work best when the 
current QMC ratio n 5 / 2 :n 7/2 was retained during the 
cheap iterations which were then used only to fix the 
total / charge rif. Also note this approach requires the 
strict equivalence of the rij calculation in both the lattice 
and impurity problems, which dictates using 

Gfo-M - |g^£^{G k M}, (12) 

in conjunction with Eq. JSJ in solving Eq. (J3J) for the bath 
Green function Gj8jj'5 mm ,' ■ 

The output of each QMC cycle is the imaginary time 
Green function G™ p (r) which contains statistical un- 
certainties, and must be Fourier transformed in order 
to define the new self-energy Hj(iuj). Following earlier 
work(2ii4& our approach is to use multipole fits 

^W-FM^wMt), (13) 

i 

using basis functions /i(r) = — e~ EiT /(e _ ^ £< + 1), which 
have Fourier transforms fi(iui) = l/(iu> — £i), thus triv- 
ially giving the overall Fourier transforms. These fits 
were carried out requiring positive weights Wji > 0, and 
constrained to give the QMC values at r = + 

F J (0+) = G; np (0+) (14) 

as well as a number of correct moments m = 1, 2, • ■ • 

p( m - 1 \o+) + F^ip-) = (-l) m Gj m) . (15) 

Here 0~ = ,8-0+ F^(t) = {d m / d,T m )F {t) , and the 
high-frequency moments are defined by G™ p (ic<j) = 

Gj^/{iw) m . In our previous work on Ce we imposed 
two (m = l, 2) of the constraints in Eq. (|15l) . however, we 
find in the present effort for Pr and Nd at the largest 
volumes that it was necessary to add a third moment 
(m = 3). The reason is that for such strongly localized 
functions where the G™ p (t) fall rapidly away from their 
t = + and j3~ values, one may get large and unphysical 
values of and G^ of opposite signs. This can lead to 
unphysical structure in the spectra far from the chemical 
potential, although we could detect no impact on the 
total energy or rij from this problem. We found that 
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adding the m = 3 or Gj constraint was sufficient to fix 
this difficulty. 

The moments Gj may be obtained from k-averages 
and appropriate traces of powers of ^lda(^) with in ad- 
dition £< 0) from Eq. © for Gf\ and both sj 0) and 

for G^ . sj 1 ' may be obtained from the double occupa- 
tion matrix^i or approximated using a Hubbard-I type 
of self-energy tuned to give the QMC rij values. 

The results reported in this paper are for a temper- 
ature of 632 K (0.004 Ry), which given the T~ 3 cost 
of the QMC calculations is about as low as is practi- 
cable. Here the imaginary time interval (0, j3) for the 
QMC was discretized by L = 80, 112, and 160 divisions, 
carrying out 10,000, 6,000, and 2,000 QMC sweeps, re- 
spectively, for each DMFT(QMC) iteration. The smaller 
number of sweeps for larger L was dictated purely by the 
L 3 expense of the calculations, and the statistical un- 
certainties necessarily increased. Nonetheless systematic 
and generally L~ 2 behavior was seen allowing extrap- 
olation to eliminate the Trotter errors^ Starting self- 
energies were already quite good, taken from converged 
DMFT calculations with the Hubbard I self-energy^ or 
from nearby volumes, or temperatures (0.01 Ry). At 
least thirty iterations were performed at each volume, 
following five discarded warmup iterations. At a num- 
ber of volumes this process was repeated taking the fi- 
nal self-energy as the initial guess and starting the pro- 
cess from the beginning. The results were unchanged 
to well within the statistical uncertainties. All Mat- 
subara sums for an argument i 7 '(iw„)e 1 "" were car- 
ried out using an asymptotic two-pole approximation 
F2 P oi(iw) — w\/(iuj — ex) + W2/{iu — £2) with the pa- 
rameters fixed by the first four high-frequency moments 
of F(iu>). The infinite sum was then given by a finite sum 
over the difference F(iuj n ) — F2 po i(iuJ n ) plus the analytic 
result for the infinite sum over F^poi, with 256 positive 
Matsubara frequencies used in the former. 

We found a noticeable increase in the QMC statisti- 
cal uncertainties for the total energy in going from Ce to 
Pr to Nd, likely reflecting the larger role played by the 
potential energy Ufd and the QMC determined value of 
the double occupation, Eq. ©, with the increasing num- 
ber of 4/ electrons. At large volumes, for example, d 
should be approximately 0, 1, and 3 for Ce, Pr, and Nd, 
respectively, at low temperatures. Similarly we found dif- 
ferences in the nature of the Trotter corrections, where 
again at larger volumes we found the coefficient of the 
L~ 2 dependence to be quite small for Ce, larger for Pr, 
and significantly larger for Nd. 

We were able to get decent spectra from the multipole 
fits Eq. l|13f) by increasing the number of poles by well 
over two orders of magnitude. As described earlier^ we 
took equally spaced £j grids of L/4 points, and systemat- 
ically eliminated poles with negative weights, examining 
O(10 4 ) such grids of varying centroid and width to find 
the best fit. For the spectra we combined the 30 best 
fits so long as the worst of these had a root-mean-square 



agreement with the QMC data no more than 20% larger 
than that of the best. More important, we also aver- 
aged over the fits for the last half of the DMFT(QMC) 
iterations. The resultant collections of O(10 3 ) poles were 
broadened by Gaussians of 0.5 eV full- width at half max- 
imum. We found systematic evolution of these spectra 
with volume providing one measure of their validity. 

C. Double occupation and moment 

Since we approximate the /-/ self-energy by the form 
Y,j mt ji m >(iu>) = Yij(iuj) Sjj'S mm >, a comparable treatment 
of the double occupation matrix is 

! n j/( 2 3 + 1 ) if j=j',m=m' 
cJ ii/[i( 2 J + 1 )] ^ j=.f, m 7^ m ' , (16) 
^5/2,7/2/48 ifjW 

noting that n 2 m =hj m , and j = 5/2 and 7/2. We obtain 
rij and djj' from the QMC auxiliary impurity problem 
by summing over the appropriate blocks in the 14 x 14 
matrix (hj m hj> m >). As these are block totals, the total 
number of / electrons per site is n/ = n 5 /2 + n 7 /2, an d 
the double occupation d appearing in the expression for 
the total energy is d = d 5 /2, 5 /2 + ^5/2, 7/2 + ^7/2, 7/2- We 
find the present djj for j = j' to be bounded above by 
their uncorrelated values from (hj m fij m i) ~ {nj m ){hj m >) 
for m 7^ m', or 

</;"' ./»; (2j • 1 ) • (17) 

Similarly, we find the djj to be bounded below by a typ- 
ical strongly correlated expression 

<% e = l j [n j -(l j +l)/2], (18) 

where lj is an integer such that lj < rij < lj+1 and Eq. 118|) 
is a piecewise linear function of rij with values rij {rij— 1)/2 
at integer rij. 

The expectation of the onsite squared / moment is 
given by (J 2 ) = 3(J|) since our Hamiltonian is rotation- 
ally invariant, and thence using Eq. i|16|) by 

(J 2 )= E V+VUih-dji)- (19) 

3=5/2,7/2 

In the local-moment regime at large volume and low 
temperature where n 5 / 2 = 1)2, and 3 for Ce, Pr, and 
Nd, respectively, n 7 / 2 — 0, and djj = rijirij — l)/2, 
Eq. jpj yields J = 5/2, 3.2749, and 7/2, respectively, 
via (J 2 ) = J(J+1). This compares to the Hund's rules 
ground state values of Jt ruc = 5/2, 4, and 9/2, respec- 
tively, and reflects the fact that including the spin-orbit 
interaction in the absence of intraatomic exchange will 
give the correct moment only when the j = 5/2 subshell 
has a single electron or hole, or is trivially empty or full. 
From the thermodynamic perspective, it is the degen- 
eracy of the Hund's rules multiplet that matters, and 
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FIG. 1: Squared local 4/ moment (J 2 ) for Ce, Pr, and Nd. 
The DMFT(QMC) results are given by the solid curves with 
data points; local-moment values, by the horizontal long- 
dashed lines; and an uncorrelated estimate with quenched 
spin orbit, by dash-dot lines. The vertical short-dashed lines 
denote transitions with significant volume collapse by shad- 
ing. Both a and 7 Ce phases are fee. 



this changes from 14!/[n!(14 — n)!] = 14,91, and 364 to 
6!/[n!(6-n)!] = 6,15, and 20 for Ce, Pr, and Nd, re- 
spectively, on adding spin orbit, which is quite an im- 
provement given the correct degeneracy 2J true + l = 6,9, 
and 10, respectively. Although not described here, we do 
find DMFT calculations with the Hubbard-I self-energy 
to yield the expected low-temperature entropy plateaus 
of In 6, In 15, and In 20 in units of fee for Ce, Pr, and Nd, 
respectively, at large volume, similar to earlier work for 
Ce without the spin-orbit interaction!^*^ 

One may use Eqs. (|17|l and (|18|l to find lower and upper 
bounds on (J 2 ) in Eq. i|19|) . If one assumes that the spin 
orbit is also quenched, rij = (2j'+l)n//14, in addition to 
the uncorrelated Eq. (|17fl . then 



(J 2 )unc =51 71/(1-^/14)74. 



(20) 



With this additional assumption, Eq. I|20|) no longer pro- 
vides a lower bound on (J 2 ) except at the smallest vol- 
umes where the spin orbit is quenched, but will never- 
theless prove useful. 



III. RESULTS 

This section reports calculations for compressed Ce, 
Pr, and Nd obtained by the LDA+DMFT method with 
a QMC determination of the self-energy, to be referred to 
more simply as just DMFT. All theoretical results have 
been carried out for an assumed fee structure and at a 
temperature of 632 K (0.004 Ry) unless otherwise noted. 
This temperature is about as low as is practical given 
the T~ 3 expense of the QMC. Nevertheless, it is cold 
enough, since previous work for Ce found that the total 
energy and entropy at 632 K are relatively close to the 
low-temperature limiti^*^ Thus, e.g., the slope of the 
7-a phase line in Ce comes primarily from the explicit T 
in the TS term of the free energy F = E — TS, and not 
from T-dependence in either the energy or entropyS 



A. Moment and spectra 

The extended aspect of the evolution from localized 
to itinerant character in the three lanthanides is illus- 
trated by the DMFT results (data points with error bars) 
in Fig. ^ for the square of the onsite 4/ electron mo- 
ment (J 2 ) calculated using Eq. (JT^J. The vertical dashed 
lines mark volumes at which transitions are observed ex- 
perimentally. The horizontal long-dashed lines at large 
volume give the strongly-correlated local-moment values 
J( J+l) with J the proper Hund's rules 5/2 value for Ce, 
although somewhat smaller than the proper values for 
Pr and Nd due to our omission of intraatomic exchange 
(see Sec. Ill C|) . It is apparent for each material that the 
local-moment regime persists under compression up to 
and through most of the stability field of the fee phase 
(large volume 7 fee phase for Ce). The dash-dot curves 
are uncorrelated approximations to (J 2 ) from Eq. I|20|) 
which also presume that the spin orbit is quenched, i.e., 
that ^5/2/^7/2 ~ 6/8, which is largely responsible for the 
significant offset of these curves from the local-moment 
lines at large volume. The DMFT results approach these 
uncorrelated and spin-orbit quenched curves under com- 
pression, reaching reasonably close agreement by the a 
phase of Ce and a-U phases of Pr and Nd. As will be 
discussed subsequently, the significance of the quenched 
spin-orbit lies in the fact that the lower Hubbard band 
is of predominant j = 5/2 character, while the 4/ Kondo 
resonance which grows at the Fermi level at the expense 
of the Hubbard side bands is of mixed j — 5/2, 7/2 char- 
acter. Therefore quenched spin orbit reflects dominance 
of the Kondo peak which may be viewed as a signature 
of itineracy. 

A quantitative measure of the degree of correlation is 
provided by the lower (strongly correlated) and upper 
(uncorrelated) bounds on the double occupation given 
by Eqs. (|17fl and (|18fl . respectively, which in turn pro- 
vide bounds on (J 2 ). The resultant curves have roughly 
the same shape as the DMFT results in Fig. 2] and 
bracket these results in each case, with the DMFT curves 
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FIG. 2: Number of 4/ electrons n s less an integer n = 1, 2, FIG. 3: Fraction of 4/ electrons of j = 5/2 character as a 
and 3 for Ce, Pr, and Nd, respectively. function of volume. 



switching from more or less perfect agreement with the 
strongly-correlated limit at large volume to much closer 
to the uncorrelated limit at small volume. Even so, we 
find in the vicinity of V ~ 12A 3 /atom that for each of 
these materials the DMFT results have switched only 
about 60% of the way from the strongly correlated to 
the totally uncorrelated limit, consistent with the fact 
that we still find residual Hubbard side bands at this vol- 
ume. On the other hand, standard paramagnetic LDA 
certainly does well enough in predicting the c/a ratio 
in bet Ce and structural characteristics of bet and a-U 
phases among the early actinides^S so it may well be 
that for all practical purposes such phases are weakly 
correlated enough. 

The moment (J 2 ) in Fig.^is a bare quantity and does 
not reflect any screening effects by the other electrons. Its 
large values at the smaller volumes are due via Eq. i|19|) 
to an increase in the number of Af electrons nf under 
compression as seen in Fig. [21 where rif — n is plotted 
with n = 1, 2, and 3 for Ce, Pr, and Nd, respectively. 
It is well known that the lanthanides undergo electronic 
s-d transition under pressure, during which the 6s and 
6p states begin to pass above the Fermi level thereby 
increasing the 5d and 4/ occupations M- Figure shows 
the fraction n 5 / 2 / n f °f 4/ electrons which are of j = 5/2 
character, and it is quite clear that there is quenching of 
the spin orbit in the region of interest in this work. 

The interplay between spin orbit and compression in- 
duced changes in 4/ spectra may be seen in Fig. where 



the total (Brillouin zone summed) but j-resolved 4/ spec- 
tra or density of states Dj (e) is plotted for Pr at a number 
of volumes. These DMFT results were obtained by multi- 
pole fits to the QMC Gj(r) as described in Sec. Ill Bl and 
were broadened by a Gaussian of 0.5 eV full width at half 
maximum. At the largest volume V = 46.0 A 3 /atom, 
one sees the pure j = 5/2 lower Hubbard band near 
—5 eV, and the mixed- j upper Hubbard band near 2 
eV, with the splitting consistent with Uf — 6.4 eV at 
this volume^ Under compression both Hubbard bands 
lose spectral weight at the expense of the growing Kondo 
peak near the chemical potential [i, which becomes dom- 
inant at the smallest volumes. Since this Fermi-level 
structure is also of mixed- j character, the population of 
these states at the expense of the initially pure j = 5/2 
lower Hubbard band relates the quenching of spin orbit 
in these materials to their growing itinerant character 
under compression. Note, however, as is especially ev- 
ident at V = 24.4 A 3 /atom, that the j = 5/2 and 7/2 
contributions to the Kondo peak are split by about the 
spin-orbit energy^ with the former more centered at fi. 
Consequently, the occupied part of the Kondo peak is ini- 
tially mostly j — 5/2 character, although the j = 7/2 part 
catches up as is also evident from Fig. [21 The Ce and 
Nd spectra are quite similar with the primary visual dif- 
ference the increasingly more prominent lower Hubbard 
band in the progression from Ce to Pr and Nd, as is to 
be expected. The Ce spectra are available elsewhere i^^^ 
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A cross check on the spectra in Fig. 0] is provided by 

D j ^)r,(p/ 7r )G j (p/2), (21) 

where this expression becomes exact in the low- 
temperature limits Equation (|21|l tracks the values at 
e = /.i in Fig. 0| to within better than 20% at all volumes. 
The volume dependence of D 5 / 2 (fJ,) from Eq. 1|21J) looks 
very much like that plotted for Ce in the case without 
spin orbit in Fig. 3 of Ref. [2J. For decreasing volume, 
these values are first quite small then rise to a maximum 
signifying growth of the Kondo peak, and then begin to 
decrease as the hybridization induced broadening of the 
4/ bands begins to dominate. Ce, Pr, and Nd all behave 
similarly with the onset of growth in the Kondo peak 
starting near the low-volume side of the fee (7-Ce) phases 
and then rising to successively higher maxima occurring 
at smaller volumes, roughly at 29, 23, and 20 A 3 /atom, 
respectively. Comparing these volumes with the tran- 
sitions indicated in Fig. QJ it may be seen that the Ce 
value reaffirms earlier observations that the initial rapid 
rise of the Kondo peak in that material coincides with 
the volume collapsei 25 i 26 ' 27 i 34 ' 35 For Pr and Nd, on the 
other hand, the growth of the Kondo peak coincides with 
the observed stability field of the dfee phase, suggesting 
a transitional role for this phase rather than being the 
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FIG. 5: (Color online) (a) Correlation energy for Pr as a 
function of volume, (b) The more interesting part of the Ce, 
Pr, and Nd correlation energies, namely the DMFT energies 
in each case less fully polarized Hartrec Fock (HF). The insets 
(offset vertical axes) show the large volume behavior more 
clearly. The ni oc in -Ehf(«1oc) refers to the number of 4/ 
bands split off below the Fermi level in the HF calculations, 
with ni oc = the customary paramagnetic solution, n\ oc > 
indicating spin- and orbitally-polarized solutions, and n\ oc — 1, 
2, and 3 for Ce, Pr, and Nd, respectively, being the fully 
polarized solutions. 



end member of the localized trivalent lanthanide series 
as has been assumed. In contrast to the j — 5/2 function, 
D 7 / 2 (n) starts to increase for decreasing volume at about 
the same place as D 5 / 2 (Ai), but does so more gradually, 
steadily growing until it becomes comparable to D 5 / 2 {n) 
at the smaller volumes. 



B. Correlation energy 

Figure OJa) shows the Pr correlation energy, namely 
the total energy less the result £hf (n\oc = 0) for a para- 
magnetic Hartree Fock solution of the effective Hamilto- 
nian Eq. jT|> . One may obtain a variety of metastable 
Hartree Fock (HF) solutions in which different numbers 
nioc of "localized" 4/ bands are split off below the Fermi 
level, ranging from the paramagnetic solution ni oc = 
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to the fully spin- and orbitally-polarized solution (maxi- 
mum 7ii oc ), which in the case of 4/ 2 Pr is ni oc = 2. The 7 
phase of 4/ 1 Ce has been described by analogous fully- 
polarized LDA+U solutions (n\ oc = l)^2iSi and the 7-a 
volume collapse by the n\ oc = 1 — > transition. 22 Fig.E^a) 
suggests a two step process in these HF solutions for Pr, 
as hrst one split-off band pops back up to the Fermi level 
under compression (the n\ oc = 2 and 1 curves cross) , and 
then later the second does the same (the rii oc = 1 and 
curves cross) , a n\ oc = 2 — > 1 — > scenario for which there 
is no experimental evidence. 

On the other hand there is some truth here, as the in- 
termediate nioc = 1 solution is really just a crude attempt 
to add 4/ spectral weight at the Fermi level in the con- 
tinued presence of some Hubbard splitting. For a truly 
correlated calculation, in contrast, there is an entirely 
continuous transfer of spectral weight from the Hubbard 
sidebands to the Kondo resonance at the Fermi level, 
consistent with the smooth behavior seen in the DMFT 
curve in Fig. Efa) . It agrees with the fully polarized HF 
result at large volumes (where the HF correctly captures 
the Hubbard splitting), but then bends smoothly away 
from this HF solution under compression (where the HF 
fails to describe the Kondo resonance) . The same behav- 
ior has been seen in earlier DMFT results for Cej^iS as 
well as in more rigorous QMC solutions for the Anderson 
lattice Hamiltoniani£2iSi It seems intuitively clear that 
such deviation from the fully polarized HF solution as 
volume is reduced should be associated with the growth 
of 4/ spectral weight at the Fermi level, i.e., the Kondo 
resonance. 

The total energy differences between DMFT and fully 
polarized HF solutions of Eq. are shown in Fig. EJb) 
for fee Ce, Pr, and Nd. Ferromagnetic order was assumed 
for the HF. The three insets (shifted vertical axes) give 
a clearer view of the large volume behavior, and the sys- 
tematically larger QMC uncertainties from Ce to Pr to 
Nd have been discussed in Sec. Ill Bl Our calculations 
of the 4/ spectra indicate growth of the Kondo reso- 
nance for each material first begins on the small-volume 
side of the experimentally observed fee or 7-Ce stability 
range (see FigOJ, and then D 5 / 2 ([i) reaches a maximum 
at ~ 29, 23, and 20A 3 /atom for Ce, Pr, and Nd, re- 
spectively. These volumes are close to where the curves 
in Fig. Efb) begin to bend away from the fully polarized 
HF (0 baseline), consistent with the association of this 
energy with the Kondo resonance. 

The fact that these deviations occur at successively 
smaller volumes from Ce to Pr to Nd is consistent with 
the ratio of the Coulomb interaction to band width 
Uf/Wf for the 4/ states, which near equilibrium volume 
(~ 34A 3 /atom) is 3.7, 4.4, and 4.9 for Ce, Pr, and Nd, 
respectively 12 Thus Pr and Nd are successively more lo- 
calized than Ce at ambient conditions, and greater com- 
pression is needed to bring about similar changes in elec- 
tron correlation. A region of negative bulk modulus is 
of course essential to obtain a first-order isostructural 
transition as in Ce. If more generally a region of low 



bulk modulus were to favor larger volume changes in 
structural phase transitions, then the delayed response 
of the softening effects in Fig.[5fb) for Pr and Nd means 
competing against the ever more dominant positive bulk 
moduli of the remaining contributions to the total energy 
as they continue to grow with decreasing volume. This 
is at least intuitively consistent with a 15% volume col- 
lapse in Ce, 9% in Pr, and apparently none in Nd. It is 
difficult to quantify this point as the QMC uncertainties 
in Fig. GJb) don't admit very reliable calculation of the 
curvature of these results. 

There are also some important differences between 
the Ce and Pr volume collapse transitions, beyond the 
fact that the former is isostructural and the latter in- 
volves a structural change. As already noted, the experi- 
mentally observed two-phase region (27.8-33.1 A 3 /atom) 
for the 7-a collapse in Ce coincides with the region of 
rapid grown in the Kondo resonance. It also overlaps 
a rather sharp breakaway of the DMFT curve from the 
HF result, as seen near 29A 3 /atom in Fig. |^b). which 
is consistent with a large negative contribution to the 
bulk modulus. Here, as for earlier results on Ce with- 
out the spin orbit as well as analyses of the An- 
derson lattice Hamiltonian, 59,60 these features corrobo- 
rate the Kondo-volume collapse model of the Ce 7-a 
transitioni^i2Si2i For Pr, on the other hand, the volume 
collapse appears to occur at pressures above the region 
of rapid growth of the Kondo resonance. The maximum 
in the Pr j= 5/2 spectral weight at the Fermi level is 
near 23 A 3 /atom (see also Fig. 0}, and the Pr curve in 
Fig-Elk) begins to bend away from the fully polarized HF 
at slightly larger volume, both of which lie outside and 
to the large- volume side of the observed two-phase region 
for Pr (19.8-21.8 A 3 /atom). If the Ce collapse serves to 
initiate the evolution from localized to itinerant in this 
material, such evolution is already well underway by the 
Pr collapse, a point which is also quite evident from be- 
havior of the local moment in Fig. ^ 



C. Equation of state 

Comparison of the theoretical and experimental en- 
ergies are given for Ce, Pr, and Nd in Figs. El [7| and 
El respectively. The HF and DMFT total energies are 
all at 632 K and for an assumed fee structure as noted. 
Experimental pressure data P(V,T) exists for Pr over a 
wide range of volume at 300 < T < 725 Kjl from which 
one may obtain the free energy F(V, T) and total en- 
ergy E(V, T) to within arbitrary constants, which we fix 
by setting E(V , T) = F(V ,T) = at the experimental 
300-K equilibrium volume Vq. Raising the temperature 
generally makes the free energy a steeper function of vol- 
ume. However, since dP/dT = dS/dV > 0, the total 
energy is generally less steep than the free energy. It is 
a fortunate coincidence that these effects roughly cancel 
so that the 300-K free energy is within 0.07 eV/atom of 
the 632-K total energy for Pr throughout the range plot- 
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FIG. 6: Total energy of Ce as a function of volume. The 
theoretical HF and DMFT results are at 632 K. The experi- 
mental energy (632 K) and free energy (299 K) are obtained 
from Refs. y and |g3 as described in the text. The theoret- 
ical bet estimate took a T — local-density bet-fee energy 
difference^ and added this to the present fee DMFT result. 
The experimentally observed 300-K equilibrium volume Vb 
and stability fields for various phases are marked. 



FIG. 7: Total energy of Pr as a function of volume. The 
theoretical HF and DMFT results are at 632 K. The experi- 
mental energy (632 K) and free energy (299 K) are obtained 
from Ref. as described in the text. The theoretical a-U es- 
timate took a T = local-density a-U— fee energy difference, 18 
and added this to the present fee DMFT result. The experi- 
mentally observed 300-K equilibrium volume Vo and stability 
fields for various phases are marked (cp denotes close packed) . 



ted for these experimental quantities in Fig. |7| noting 
again that both are zeroed at Vq. Experimental mea- 
surements of P(V, T) for Ce exist up 208 GPa at 300 K£ 
and over the temperature range 299-573 K within the 
fee (7 and a) regime^ Assuming dP/dT is independent 
of volume, we may use the former results to extend the 
latter up to higher pressures. Extrapolating these results 
also to somewhat higher temperatures, one may obtain 
the F(V, 299 K) and E(V, 632 K) curves seen in Fig. © 
which like Pr are also fairly close. We presume that it 
is also legitimate to compare room-temperature free en- 
ergies for Nd to our 632-K calculated total energies for 
that material. 

The Debye temperatures for Ce, Pr, and Nd are all be- 
low 150 K at ambient conditions^ and with reasonable 
Griineisen parameters all three materials should still be 
in the high-temperature limit at 632 K throughout the 
volume range studied here. The phonon contribution to 
the total energy is then 3fceT and has no impact on the 
volume dependence examined in Figs. HjHSJ 

The total energies in Figs. HJIEl provide another per- 
spective on the correlation issues discussed earlier. The 
fully-polarized HF results (dashed) curves are analogous 
to LDA+U and should do well at large volumes in the 



strongly localized limit as is the case here. The DMFT 
results (open circles) then show how such static mean 
field theories begin to break down as volume is reduced 
due to their inability to account for the growing Kondo 
resonance and its contribution to the correlation energy. 
This effect is most pronounced for Ce in Fig. EI less so for 
Pr in Fig. [7| and smallest for Nd in Fig. [SJ reflecting the 
fact that Pr and Nd start out successively further in the 
localized limit than Ce, therefore requiring greater com- 
pression to achieve comparable changes in the correlation 
energy. Although the agreement between the DMFT re- 
sults and the experimental isotherms is not ideal, and 
will be discussed further, it is clear that experiment does 
confirm this systematic progression. 

A comparison of theoretical and experimental values 
for the bulk properties is given in Table [Q The theory 
results were obtained from 6-term fits£& to the energies 
over the range 17-49 A 3 /atom in order to better aver- 
age over scatter in the QMC results. Table|Jshows quite 
decent agreement with experiment^ and is a partial val- 
idation of the present effective Hamiltonians Eq. JJJ and 
the manner of total energy calculation Eq. Q . Of par- 
ticular note is the fact that the DMFT results for Vq and 
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FIG. 8: Total energy of Nd as a function of volume. The theo- 
retical HF and DMFT results are at 632 K. The experimental 
free energy (300 K) was obtained from Ref. Il3l The experi- 
mentally observed 300-K equilibrium volume Vb and stability 
fields for various phases are marked (cp denotes close packed) . 



TABLE I: Comparison of theoretical and experimental val- 
ues for the equilibrium volumes Vb(A 3 /atom) and bulk mod- 
uli Bo(GPa). The theory used fits to the 632-K results in 
Figs.|SJ[H| while the experiment is at 300 K. However, the Ce 
HF results change by only 0.4% and —2.6% for Vb and Bo, 
respectively, on reducing the temperature from 632 to 316 K. 
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Bq are in better agreement with experiment and 3% and 
36% smaller, respectively, than the fully polarized HF 
values in the case of Ce, as compared to ~ 1% and 6- 
9% smaller, respectively, for Pr and Nd. This suggests a 
small but not unimportant effect of the Kondo resonance 
even at P = in Ce, effects which are shifted to smaller 
volume and have less impact for Pr and Nd. 

An estimate of the structural contribution to the en- 
ergy at smaller volumes is provided by the "bet" and 
"a-U" results (x symbols) in Figs. HO and [7| for Ce and 
Pr, respectively. These were obtained by adding para- 
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FIG. 9: (Online color) Experimental 300-K pressure-volume 
data for (a) Ce^ (b) PrjLil&L and (c) Nd t 13 i 68 ' 69 as com- 
pared to theory. The observed volume collapse transitions 
in Ce and Pr are shaded. The thin solid lines give T — 
LDA results^ while the HF, DMFT, and structural estimates 
("bet" Ce, "a-U" Pr) are negative volume derivatives of the 
total energies in Figs. |SJ|H| 



magnetic local-density functional values for the bet-fee 
energy difference in CejSi and the a-U-fcc difference in 
Prji& to the DMFT fee values. Throughout the ranges 
shown in Figs.|HlandElthese structural energy differences 
are 10% or less for Ce and 25% or less for Pr of the more 
dominant correlation energy contribution represented by 
the separation between the DMFT and the polarized HF 
curves. These structural corrections improve the agree- 
ment with experiment somewhat for Ce and Pr, however, 
a similar correction for Nd would move the DMFT results 
to lower energies further away from the experiment. 

Figure ED gives the pressure- volume results correspond- 
ing to Figs. EHH1 with experimental 300-K data shown 
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for Ce^ Pr^^ an d Nd, 13 i 68 i 69 Except for the T = 
LDA results (thin solid lines) r& the theoretical results 
are negative volume derivatives of the 632-K total ener- 
gies in these figures, an approximation justified by the 
-Fexp(300K) ~ i? CX p(632K) comparisons in Figs. and 
Multi-term fits to the energies were used in obtaining 
the derivatives^ The DMFT points (open circles) were 
fit ignoring the error bars, with separate fits made for the 
regions above and below the observed collapse transitions 
for Ce and Pr. It is evident that the present paramag- 
netic DMFT results do well in both low-pressure localized 
and high- pressure itinerant extremes, especially consid- 
ering the structural corrections [dashed lines in Figs.^a) 
and (b)] in the later regime, although there are evident 
discrepancies in between as will be discussed. Stan- 
dard paramagnetic LDA does well in the itinerant regime 
as seen by the thin black lines, but must be combined 
with one of the magnetically-ordered modified-LDA tech- 
niques (e.g., LDA+U) to capture the larger volume be- 
havior. This then leads to a magnetic order-disorder 
transition which is not observed at room temperature, 
and, e.g., an incorrect prediction of a volume collapse 
(16% at 39 GPa) in Nd^ 

Turning to the collapse transitions, the kink in the 
DMFT total energy for Ce near 29 A 3 /atom in Fig. © 
lies within the experimental two-phase region, is within 
the QMC error bars, and is similar to earlier results with- 
out spin orbiti^' "' 3 ^' The addition of the spin-orbit inter- 
action has lowered the 7-phase energy more than that 
of the a phase, so that a common tangent construction 
to the DMFT energy curve in Fig. gives a pressure of 
4 ± 2 GPa as compared to a slightly negative value in the 
earlier work. While the experimental transition ranges 
from (an extrapolated) —0.7 GPa at T — to 1.8 GPa at 
485 K, 62 an error of 4 GPa is not huge on the scale of 
Fig.EOa). The volume change at 4 GPa in the DMFT re- 
sults (thick solid lines) in Fig. |3f a) is in decent agreement 
with that observed at room temperature. In the case of 
Pr, the estimated a-U curve (dashed line) in Fig. Etb) is 
in good agreement with experiment over all but the low- 
est observed range of this phase. This raises the question 
of whether the structural energy difference may be criti- 
cal to the Pr collapse. Indeed, while we find a softening 
in the Pr equation of state for the fee structure due to the 
contribution in Fig.|SJb), there is no convincing evidence 
that this drives the bulk modulus negative in the vicinity 
of the collapse. 

Finally, we turn to the discrepancies between the- 
ory and experiment in Figs. EHH1 The fact that the 
DMFT results for Ce in Fig. are too high at smaller 
volumes in comparison with experiment, those for Pr 
in Fig. are reasonable, while those for Nd in Fig. |H1 
are somewhat low is suggestive of the intraatomic ex- 
change interaction, K, which has been omitted in the 
present work. In the strongly localized regime includ- 
ing V ~ Vo, this exchange interaction should shift the 
energy by 0, —K, and — 3K for Ce, Pr, and Nd, re- 
spectively. At small volumes, the itinerant expectation 



would be ~K J2<rI2m<m'( n ™<r n ™'<r) ~ -3Knj/U tak- 
ing (n mrT n m ' a) ~ (nf/14) 2 . If the total energies are 
shifted to at Vb as done in these figures, this would sug- 
gest shifts of -0.31, 0.08, and 0.69 eV at V=15 A 3 /atom 
for Ce, Pr, and Nd, respectively, taking values of n/ from 
Fig-Eland K ~ 1 eV. This would imply completely itiner- 
ant states at 15 A 3 /atom which is unlikely especially for 
Nd. Indeed, polarized HF calculations suggest exchange 
corrections at this volume of —0.13, 0.07, and 0.25 eV, 71 
respectively, although this also is only an estimate. Nev- 
ertheless, these estimates do suggest that intraatomic ex- 
change would significantly improve the present compar- 
isons between theory and experiment. Unfortunately, a 
rigorous inclusion of the exchange interaction, one that 
would also yield the the correct Hund's rules values for 
the Pr and Nd local moments, would require DMFT cal- 
culations for the full /-/ Coulomb interaction including 
all four Slater integrals and the non-density-density ex- 
change and pair hopping terms. This has been a chal- 
lenge for QMC, although there is recent progress^ Sim- 
ilarly, implementing mutual self-consistency between the 
LDA and DMFT parts of the calculations, as has been 
advocated^ would be prohibitively expense with the 
present QMC implementation of the self-energy. Both 
improvements as well as DMFT calculations of the struc- 
tural energy differences could be considered with faster 
although less rigorous approaches to the self energyiS 

IV. SUMMARY AND DISCUSSION 

A better understanding of the electron-correlation 
driven volume collapse transitions in the compressed lan- 
thanides may come from putting this behavior in context 
of their extended evolution from localized to itinerant 
character. To this end, the present paper has reported 
calculations for compressed Ce (4/ 1 ), Pr (4/ 2 ), and Nd 
(4/ 3 ) using a combination of the local density approxima- 
tion (LDA) and dynamical mean field theory (DMFT), 
so called LDA+DMFT. Results for the 4/ moment (J 2 ), 
spectra, correlation energy, and equation of state among 
other quantities have been presented over a wide range 
of volume at a temperature of 632 K. This temperature 
is the lowest feasible with the present quantum Monte 
Carlo (QMC) implementation of the self energy, yet is 
still reasonably close to the low-temperature limit 
While a face-centered cubic (fee) structure was assumed, 
LDA estimates of the important structural energy dif- 
ferences are significantly smaller than the relevant con- 
tributions to the correlation energy which may then be 
described as leading order volume-dependent effects. 

We find the three lanthanides to remain rather strongly 
localized under compression from ambient conditions up 
through the observed stability fields of the fee (7 Ce) 
phases, in the sense that the 4/ moments are close to 
the Hund's rules values, there are fully formed Hub- 
bard sidebands which are themselves a signature of the 
local moments, and little 4/ spectral weight lying in 
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between at the Fermi level. Subsequent compression 
brings about significant deviation of the moments from 
the Hund's rules values, growth of 4/ spectral weight 
at the Fermi level (the Kondo resonance) at the ex- 
pense of the Hubbard sidebands, an associated soften- 
ing in the total energy, and quenching of the spin orbit 
given that the Kondo peak is of mixed- j character in 
contrast to the predominantly j — 5/2 lower Hubbard 
band. The most dramatic evolution in these signatures 
is seen to coincide with the two-phase region of the 7— a 
phase transition in the case of Ce, consistent with earlier 
wor ] t! 33 i 34 i 35 an( j m agreement with the Kondo volume- 
collapse scenario For Pr and Nd, on the other 
hand, these signatures change most rapidly over the vol- 
ume range where the distorted fee (dfee) structure is ex- 
perimentally observed to be stable, suggesting that this 
phase is transitional and not part of the localized triva- 
lent lanthanide sequence. 

Only on subsequent compression is Pr experimentally 
observed to undergo a collapse from the dfee phase to an 
a-U structure, while Nd passes from dfee to a-U through 
two other low symmetry phases without any substantial 
volume changes. Due to the increasing but incompletely 
screened nuclear charge, the lanthanides shift towards 
the localized limit for larger atomic number, and we see 
a similar off-set to smaller volume of the above mentioned 
signatures from Ce to Pr, and then Nd. In particular, a 
softening contribution to the total energy associated with 
growth of the Kondo resonance must compete with the 
remaining contributions which become ever more domi- 
nant due to a steadily increasing bulk modulus as vol- 
ume is reduced. If a region of low bulk modulus were 
to favor larger volume changes in structural transitions, 
then this would be qualitatively consistent with the ob- 
served sequence 15% (Ce), 9% (Pr), and none (Nd). This 
speculation is apparently contradicted by the 5% and 6% 
collapse transitions in Gd (4/ 7 ) and Dy (4/ 9 ), respec- 
tively. However, these lanthanides correspond to filling 
a different spin-orbit subshell, which may significantly 
complicate matters given the profound manner in which 
spin-orbit is involved. 

There are still some correlation effects evident at even 
the smallest volumes considered here, such as residual 
Hubbard side bands in the 4/ spectra and moments 
whose values have evolved only ~ 60% of the way from 



the strongly localized Hund's rules values to those char- 
acteristic of totally uncorrelated electrons. Nevertheless, 
standard paramagnetic LDA does quite well within the 
bet phase of Ce and the a-U phases of Pr and Nd, and 
so for all practical purposes these phases appear weakly 
correlated enough. The success of LDA calculations for 
the structural parameters in these phases for Ce as well 
as the light actinides is also well known^Si 

This work has also included a detailed quantitative 
comparison between the present LDA+DMFT results 
and experiment for the total energies and pressures of 
Ce, Pr, and Nd. The comparison is encouraging, and 
serves to corroborate such theoretical observations as the 
systematic offset of Kondo-like correlation signatures to 
smaller volume for increasing atomic number in these 
three lanthanides. There are also clear deficiencies, most 
notably the need to include intraatomic exchange, ide- 
ally in its full rotationally invariant form to enable rig- 
orous calculation of the Hund's rules moments in general 
cases. The need for accurate LDA+DMFT structural 
energy differences is also apparent. These improvements 
as well as making the LDA and DMFT parts mutually 
self-consistent will require a self-energy treatment that is 
both more precise and considerably less expensive, which 
unfortunately is likely to require giving up some of the 
rigor of the present QMC approach. 
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